One-dimensional Ising model with long-range and random 

short-range interactions 

A. P. Vieira^ and L. L. Gongalves* 

Departamento de Fisica da UFC, Cx. Postal 6030 
60451-970 Fortaleza (CE), Brazil 

t Present address: Institute- de Ffsica da USP, 05315-970 Sao Paulo (SP), Brazil. 

' E-mail: apvieira@if.usp.br 
* E-mail: lindberg@fisica.ufc.br 

February 1, 2008 



Abstract 

The one-dimensional Ising model in an external magnetic field with uniform long-range interactions 
and random short-range interactions satisfying bimodal annealed distributions is studied. This general- 
izes the random model discussed by Paladin et al. (J. Phys. I France 4, 1994, p. 1597). Exact results 
are obtained for the thermodynamic functions at arbitrary temperatures, and special attention is given 
to the induced and spontaneous magnetization. At low temperatures the system can exist in a "ferri- 
magnetic" phase with magnetization < a < 1, in addition to the usual paramagnetic, ferromagnetic 
and antiferromagnetic phases. For a fixed distribution of the random variables the system presents up to 
three tricritical points for different intensities of the long-range interactions. Field-temperature diagrams 
can present up to four critical points. 



1 Introduction 

The one-dimensional Ising model with nearest-neighbor and uniform infinite-range interactions was intro- 
duced by Nagle Q| and presents interesting behaviour such as the existence of first- and second-order para- 
magnetic transitions. The crossover between short- and long-range interactions in the model was later re- 
examined by Kislinsky and Yukalov ^| , and the study of the model for arbitrary interactions was presented 
by Vieira and Gongalves || . The model is equivalent to a linear chain approximation of higher-dimensional 
models, and it has been successfully applied to the study of quasi-one-dimensional materials |5|. 

The presence of infinite-range interactions is responsible for the existence of ferromagnetic ordering at 
finite temperatures, and is expected to induce even richer behaviour in the presence of disorder. This 
is confirmed by the analysis of the quenched site-dilute system by Slotte ||, who found that the system 
exhibits 11 topologically different field-temperature phase diagrams as p and a are varied, where p is the 
dilution concentration and a measures the ratio between the intensities of nearest-neighbor and infinite- 
range interactions. Recently, Paladin et al. studied a random version of the model, at zero field, in 
which nearest-neighbor interactions can assume opposite values ± J with equal probability. In the quenched 
limit they showed that, as a is varied, infinitely many "ferrimagnetic" ground state phases, with average 
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spontaneous magnetization < itq < 1, are observed. Because of the difficulties involved in the mathematical 
treatment of the system at arbitrary temperatures, they also considered the annealed limit, which presents 
simpler ground state behaviour, but where a ferrimagnetic phase is still present. 

In a recent paper || the present authors considered a generalization of the annealed model studied 
by Paladin et al. by introducing an external field and random nearest-neighbor interactions which satisfy 
bimodal annealed distributions. The exact solution of the model was presented, and results for the bond- 
dilute case were briefly discussed. This paper considerably extends the results of the previous one by 
considering all possible cases for the distribution of the random variables. The exact solution of the model 
is reviewed in Sec. || by using the approach of Thorpe and Beeman ^ , which maps the system onto a 
regular Ising model with an effective temperature dependent interaction. The ground state properties of 
the system are discussed in Sec. |[ and it is shown that ferrimagnetic phases do appear for a very large 
range of distributions. At finite temperature the system presents a very rich behaviour, with many phases, 
as discussed in Sec. ^. Finally, in Sec. ^ the main results of the paper are summarized. 



2 The model 

Consider a closed-chain (TV sites) Ising model in an external magnetic field with uniform long-range interac- 
tions and random nearest-neighbor (n.n.) interactions satisfying the bimodal annealed distribution 

p(Kj) = p5(nj - J A ) + (1 -p)S(k j - J B ). (1) 

The Hamiltonian of the model is written as 

N j N N 

H = ^2 '•V T .'' T .' • i - ]y ^2 u i a i ( 2 ) 

which generalizes both the pure model UUP and the random model studied by Paladin et al. 0, 
corresponding to h = 0, J a = —Jb = J and p = h- 

The exact solution of the model was given by the present authors in a previous paper ^|, where an 
even more general model was considered, for which long-range interactions satisfy a strongly correlated 
distribution. Alternatively, the thermodynamic functions of the model can be easily obtained from the 
general approach of Thorpe and Beeman [^| . This consists of performing the partial trace over the disorder 
variables, mapping the system onto a regular Ising model with an effective n.n. interaction K, whose 
properties are well known Q |^, [|. F° r the random system, the free energy per spin as a function of the 
average magnetization a is given by 

f a (a) = f(K, a) - k B T {pin [cosh(/f - (3J A ) - e{K) sinh(X - 0J A )] + 

+ (l-p) In [cosh(X - f3J B ) - e{K) sinh(X - (3J B )} - [p\np + (1 - p) ln(l - p)}} , (3) 
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where 



f{K, a) = -k B T\n 



e K 



cosh h + ( sinh 2 h + e~ iK 



lo\ (4) 



with 



h = /3h + 2f3Ia = h + 21a, (5) 
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is the free energy of the regular model [0 with (uniform) effective n.n. interaction K, determined from the 
solution of the equation 

P ' =0, (6) 



coth(AT - (3J A ) - e(K) coih{K - (3J B ) - e(K) 
and 

2e- 4K (smh 2 h + e~ 4K ^ < 



e(K) = (a Jtrj+1 ) K = 1 3 ^ (7) 

cosh ft + ( sinh 2 h + e- iK 



is the nearest-neighbor spin-spin correlation function of the regular model. As usual T is the absolute tem- 
perature, ks is Boltzmann's constant and /3 _1 — ksT. The equilibrium value of the average magnetization 
a is that which minimizes / (<j), whose relative extrema are determined from the solution of 



dfa 

da 



= 

h.T 



sinh h , , 



' sinh 2 h 



The equation of state can be calculated by determining, among the solutions of Eq.|§|, the absolute minimum 
of f a at given field h and temperature T. 

Competition between short- and long-range interactions is responsible in the pure model for the existence 
of first order phase transitions at both zero and non-zero field. It is then natural to expect that it will produce 
even richer behaviour in the presence of disorder. Assuming the long-range interaction to be ferromagnetic 
(I > 0) competition between interactions is obtained by considering Jg bonds as antiferromagnetic ( Jb < 0) 
and allowing J a bonds to have arbitrary character (Jb < J a < +oo). These are the cases considered in this 
work. As shown for the pure model j| , for / < the system is fully frustrated and so no spontaneous order 
is possible even at T = 0. 



3 Ground state properties 

Since at T = the internal energy E and the free energy are equal, the ground state properties of the 
model can be obtained by looking for those configurations which minimize E for fixed values of the various 
parameters. For annealed distributions the random variables can adjust themselves so as to minimize the 
contribution to the free energy due to n.n. interactions, with the sole restriction that the concentration of 
different kinds of bonds be satisfied. This is in general achieved at fixed average magnetization a by the 
formation of interaction domains fi |^| , which avoids frustration effects. On the other hand, the contributions 
to the internal energy due to long-range interactions I and to the external field h depend only on a and are 
given by —la 2 and —ha, respectively. 

The possible situations are analyzed below. 



3.1 Case J B < J a < 

This is the general case of antiferromagnetic n.n. interactions, due to the symmetry of the model with 
respect to the transformation p — > 1 — p, J a <-> Jb- When both long-range interactions I and field h have 
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small intensities compared to those of J a and Jb the system orders antiferromagnetically, satisfying all n.n. 
bonds. That structure has its energy per site in the thermodynamic limit given by 



Eaf = P-Ja + (1 - p)Jb, 



(9) 



where p is the concentration of J a bonds. For intermediate intensities of / and h, ferromagnetic ordering may 
be induced in the J a domains, while Jb bonds are still strong enough to maintain local antiferromagnetic 
ordering. It is easy to verify that such structure has an average magnetization a = p, and was called by 
Paladin et al. the "ferrimagnetic" structure JtJ , whose energy is given by 



-pJa + (1 - p)Jb - Ip 2 - hp. 



(10) 



For I or h strong enough the system orders ferromagnetically throughout. Such structure, in which no 
nearest-neighbor bond is satisfied, has unit average magnetization and energy given by 



-pJ A - (1 - p) Jb - I -h. 



(11) 



To determine the range of parameters in which the various structures correspond to equilibrium one has to 
compare the relative values of Eaf, Ep err i and Ep erro . Bearing in mind that for this kind of system a and 
h must have the same sign, one finds that the T = isotherm is given by 



J B < J a < \ Jb and / < 2{J A - Jb) = I 



0, if < \h\ < hi; 
p, if hi < \h\ < h 2 ] for < 

1, if \h\ > h 2 , 



J A > \ Jb and 



Ja = 
Jb-Ja - 



fi 



p* and I < I 



P< 
or 

p > p* and I < -\J A 



fi 



0, if < \h\ < h 3 ; 

1, if \h\ > h 3 , 

p, if < \h\ < h 2 ; 
1, if \h\ > h%, 



Jb < Ja < \ Jb or J a > \ Jb, P < P* 
for < and 

2(J A -J b )<K -2[ P J a + (1 -p)Jb] 



(12) 



for J a > h Jb, P > P* and - | J A < I < 



2 j ■ 



1+p 



1, if \h\ > 0, for < 



where 



Jb < Ja< \ Jb and / > -2[pJ A + (1 - p) Jb] 
or 

( p < p* and I > -2\pJ A + (1 - p)J B ] 
Ja > \ Jb and < or 

{ p> p* and / > - j^Jb 

hi = -2J A - pi, 
h 2 = -2J B -(l+p)I, 
h 3 = -2\pJ A + (l-p)J B ]-I. 



(13) 



3.2 Case J B < < J A 

In this case, even for Ja — and /, h —> + internal energy is minimized by a structure where antiferromag- 
netic ordering prevails in J B domains while the rest of the system orders ferromagnetically. This corresponds 
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to the ferrimagnetic phase (a = p), and so for J a > there is no stable antiferromagnetic structure. Again 
for I or h strong enough the system orders ferromagnetically throughout (a = 1). The relevant internal 
energies per spin in the thermodynamic limit are 

E F erri = ~pJ A + (1 - P)Jb ~ Ip 2 ~ hp (14) 

and 



-pJ A -(l-p)J B -I-h, (15) 



so that the T = isotherm is given by 



for I < 



2J B 



p, if < \h\ < -2J B - (l+p)I 

1, if |/i|>-2J B -(l+p)J ' 1U1 ^ (i+P)' 

(16) 



1, if \h\ > for I > - 



2Jb 



(T+P)' 

4 Finite temperature properties 

From the results of the last section, one concludes that there are essentially three different regimes for the 
ground state behaviour of the system at zero field, depending on the distribution of the random variables. 
For Jb < Ja < \Jb the spontaneous magnetization (oo = cr|/j = o) can assume only the values <jq = and 
Co = 1, corresponding to antiferromagnetic and ferromagnetic phases, respectively. For \J B < Ja < 
ferrimagnetic (<7o = p) phases become possible, in addition to antiferromagnetic and ferromagnetic phases. 
Finally, for J a > only ferrimagnetic and ferromagnetic phases exist in equilibrium. 

The ground state properties lead to the conclusion that, at finite temperatures, the free energy as a 
function of cro presents up to five relative minima in the interval — 1 < ao < 1, corresponding to the 
symmetric ferromagnetic (<ro — ±1) and ferrimagnetic (<to — ±p) phases and to the antiferromagnetic (or 
paramagnetic) phase (ao = 0). In the presence of an external field h > the relative minima of the free 
energy / (<x) are shifted towards the positive a axis and the values of /a(c) at those minima are altered, 
which may result in discontinuities or critical points in the a x h isotherms. 

The second order paramagnetic transition temperature T c is obtained, as in the pure model S , by taking 
the a — > limit with h = in EqJ|, noting however that the effective interaction K depends implicitly on a 
through e(a). Therefore, one obtains 

e- 2Ko = e- 2K \ a=Q =2f3 c I, (17) 
where = ksT c . Taking the a — * limit in Eqs.^j and |t] one also obtains 

tannic =ptanh(/3 c JA) + (1 - p) tanh(/3 c J B ), (18) 
and by combining the two previous equations one gets 

1 - ptanh(/3 c J A ) - (1 - p) tanh(/3 c J B ) 



1 + ptanh(/5 c J A ) + (1 -p) tanh(/3 c J B ) 



= 2/3 c I, (19) 



which directly determines T c . 

As discussed for the pure model jj), the tricritical point can be found by solving the system 



dfa 

da 



0, /a«)=/a(0), (20) 
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and imposing that the a t — ^ root of the second equation be quadruple by means of a Taylor expansion, 
which is equivalent to requiring that all derivatives up to fourth order of f a with respect to a be zero at 
(7 = 0. Taking into account the implicit dependence of K on a, one gets 

-4K _ 3 + QK n = 0j ( 21 ) 



e 



where 



K" 



d 2 K 



8a 2 



i2 



a=0 



8 ,pjs2 p(l - p) [tanh {(3 J A ) - tanh ((3J B )] ^ 

{l - \p tanh (/3J A ) + (1 -p)tanh(/5J B )] 2 } 

For fixed J a, Jb and p, by solving the system composed of Eqs.[l8|, |l] and |2^ the coordinates It r and T tr of 
the tricritical point are determined. In the non-random limit (p = 0, p = 1 or J a = J B ) Eq.|2^ gives Kq = 
and one obtains the pure model result [0, H || 

e-a/W = 2 /3 tr I tr , 

e -4/w _ 3 = o 

The following subsections analyze the behaviour of the system for different values of J a and Jb, using 
the renormalized parameters 

5=lA a =±_ ..- h e- kBT (24) 
Jb Jb \Jb\ \Jb\ 

4.1 Case J B < J A < \J B < (± < <5 < 1) 

For this case and at zero field (/i = 0) the behaviour of the system is qualitatively the same as for the pure 
model, regardless of the value of p, since there is no stable ferrimagnetic phase. For very strong long-range 
interactions / the system is ferromagnetic at low temperatures and there is a second order paramagnetic 
transition. As / gets weaker, passing through the tricritical point / = It r paramagnetic transitions change 
to first order and, for / < If err o = ~2 [pJa + (1 — p)Jb] the ground state becomes antiferromagnetic, and 
<j = for all T. 

In the presence of a magnetic field h > and for / < Ifi = 2(Ja — Jb) the a x h isotherms present two 
discontinuities at T — (from a — to a = p and from a — p to a = 1), which persist at low temperatures, 
inducing qualitative differences compared to the pure model. Figures 1(a) and (c) show a x 7 isotherms for 
various values of the renormalized temperature 9 = ksT/ \Jb\ when S — | and p = 0.5, with a = —0.5 and 
a = —0.6 , respectively, while corresponding 7 x 9 phase diagrams are shown in Figs. 1(b) and (d). In both 
cases there are two discontinuities in the T = isotherm for h > (or, by symmetry, h < 0), producing four 
first order lines in the phase diagrams. For a — —0.5 those lines remain separate, ending at critical points 
corresponding to temperatures 9hc and 9' hc , but for a — —0.6 the lines can intersect, giving rise to other two 
symmetric lines which terminate at a critical point. In fact, one finds that for othcr2 — —0.5751 the critical 
points corresponding to 9 hc touches the other first order lines, which produces for a < cthcri the behaviour 
observed in Figs. 1(c) and (d). 

Figure 2 shows the 9 x a diagram for 8 = | and p — 0.5. At zero field, as already stated, the behaviour 
is similar to that of the pure model, with the tricritical point separating regions of first and second order 
paramagnetic transitions. At non-zero field, on the other hand, there are two critical temperature (9h c ) 
curves for a — * 0~. One of the curves, as in the pure model, ends at the tricritical point, giving rise to the 
second order transition temperature 9 C . The other curve ends at point Phcr2, reaching the 9h cr 3 curve, which 
touches the a axis at a/i = — |. 
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4.2 Case \J B < J A < (0 < 8 < ±) 

The most interesting kind of behaviour occurs for values of S G (0, ^). The existence of ferrimagnetic phases 
at zero field for p > p* = Ja/(Jb — J a) is responsible for strong qualitative changes in the phase diagrams. 

For p > p* and / > If err i = ~2Ja/p the system can undergo first order ferri- ferromagnetic transitions 
as the temperature is raised, as observed in the x a diagrams of Figs. 3 and 4 for 6 = | and p — 0.6 and 
p = 0.7, respectively. In both figures dotted curves show first-order ferri-ferromagnetic temperatures Of. For 
p* = 0.5 < p = 0.6 < p cr 2 — 0.6596 there is coexistence of paramagnetic, ferrimagnetic and ferromagnetic 
phases at point P cr 5, where the Of curve reach the first-order paramagnetic transition temperature t . So, 
to the left of P cr 3 first-order ferri-paramagnetic transition occur at t . On the contrary, for p — 0.7 > p cr 2 
the Of curve does not reach t , ending at point P cr together with 0' hc . 

The behaviour of the spontaneous magnetization can be understood by analyzing the evolution of 
the free energy / a (co) as the temperature is varied. For 5 — ■k and p = 0.7 Figs. 5(a) and (b) show 
/(do) = /a(ffo)/ \Jb\ at different temperatures for a — —1.10 > a cr ~ —1.1018 and a = —1.13, respectively. 
In both cases the ground state is antiferromagnetic (|<7o| — p) and there exist metastable states for | ctq ] = 1 
and do = 0. For a = —1.10 the metastable ferromagnetic state disappears at ~ 0.2 and the spontaneous 
magnetization increases with temperature in a certain range. At 0t — 0.430 the system undergoes a first 
order ferri-paramagnetic transition. On the other hand, for a = —1.13 first-order ferri-ferromagnetic and 
ferro-paramagnetic transitions occurs at Of ~ 0.115 and t ~ 0.532. 

Zero field transitions can be better observed in Figs. 6 and 7, which show Co as a function of the 
renormalized temperature for the case 5 — ^ with p = 0.6 and p — 0.7, respectively. For p = 0.6 
the ground state is ferromagnetic for a < ctf erro = —1.25, and the system undergoes first-order ferro- 
paramagnetic transitions for a > a tr — —1.6997. If af erro < a < af erri ~ —1.1111 the ground state 
is ferrimagnetic, and if a < a cr ^ ~ —1.2009 both first-order ferri-ferromagnetic and ferro-paramagnetic 
transitions occur. If a > a cr ^ only first-order ferri-paramagnetic transitions are observed. For p = 0.7, 
when atr — —1.4960, af erro ~ —1.1764 and a/ err i ~ —0.9523, the behaviour is much the same, with the 
important exception that ferri-ferromagnetic transition occur if aferro < a < a cr ~ —1.1018, and exactly if 
a = a cr there is a point in the ero x curve where Co is a continuous function, but its thermal derivative 
diverges. An important feature of this model that is evident in Fig. 7 is the increase of cto with temperature 
in a certain range, which is verified in some perovskite-like materials such as Zni-xGa^MngC (x < 1) and 
CuMnaN Finally it should be noted that for \a\ greater than what is shown in the figures the behaviour 
of do becomes identical to that of the pure model, with the ferro-paramagnetic transition changing from first 
to second order at the tricritical point. 

The effects of the existence of ferrimagnetic phases at zero field on the 7x6* phase diagrams can be 
observed for S — | with p — 0.6 and p = 0.7 in Figs. 8(a)-(f) and Figs. 9(b) and (d). For p — 0.6 figures 
8(a)-(f) show the evolution of 7 x diagrams from a = — 1 to a = —1.21. If a > af err i — —1.1111 the 
diagrams are similar to those shown in Fig. 1(b). Exactly at a = af err i there appears a zero- field first-order 
line, which splits at = Of, giving rise to two symmetric first-order lines. If a > ahcr2 — —1.1295 symmetric 
lines intersect at = 0hcr3, and 7hcr3, the corresponding field, becomes zero at a cr 3 ~ —1.2009. For p = 0.7 
there is no intersection of symmetric first-order lines, as shown in Figs. 9(b) and (d). Figures 9(c) and (d) 
illustrate the behaviour of the system if a — a cr , when one of the coexistence curves just touches the a axis. 

As the value of S becomes smaller a reentrant behaviour appears in the ferri-ferromagnetic transition 
temperature {Of) lines, as can be seen in Fig. 10, which shows x a diagrams at zero field for 5 = ^ and 
various values of p. The reentrant behaviour of 6*/, evident in Figs. 10(b) and (c), is reflected in the thermal 
behaviour of the spontaneous magnetization ero shown in Fig. 11, where two first-order ferri-ferromagnetic 
transition occur for a cr ~ —1.3108 < a < a re — —1.2710. The existence of reentrant behaviour in random 
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models with competing interactions is common in both annealed Q and quenched |TlJ systems. 

Reentrant behaviour in 9 x a diagrams is related to the structure of first-order lines in 7 x 9 phase 
diagrams, as observed in Figs. 12(b) and (d), for 8 = |, p — 0.5, and a = a re ~ —1.2710 and a = —1.29, 
respectively. In the first case the symmetric upper and lower first order lines are tangent to the 7 = axis 
at 9 — 9 re , while in the second case those lines are composed of two distinct segments, touching the 7 = 
axis at temperatures 9f and 0'j, producing the discontinuities observed in the corresponding (Tq{9) curve in 
Fig. 11. 

For < 8 < 0.1649 the 9 x a diagrams can present up to three tricritical points for p S (p' rt ,p r t). This 
is exemplified for 5 — 0.16 and p' rt ~ 0.4504 < p = 0.453 < p r t — 0.4567 in Fig. 13. The existence of 
three tricritical points (Ptr, P't r ancl P'tr) * s reflected in the appearance of another zero-field second-order 
paramagnetic transition line for a' t ' r < a < a' tr . That line changes from only one point for p — p' rt , where 
P' tr and P t " coincide, until it joins the usual 9 C line for p = p rt , where P t " coincides with P tr . An example 
of 9 x a diagram with two tricritical points is given in Fig. 14. It should be noted that point P C 2, which 
marks the intersection of the 9f, 9' c and 9 t lines, is not a tricritical point, since all those lines correspond to 
h = 0, and so the critical exponents at P c2 have the usual mean-field values (a = 0, /3 = |, 7 = 1, 8 = 3, 
within the usual convention of critical indices). 

4.3 Case J A > (6 < 0) 

In this situation there is no antiferromagnetic ground state, as discussed in Sec. ||. Irrespective of the value 
of a = I / I Jb| there is always a second-order ferri- or ferro-paramagnetic transition at a finite temperature. 
This is illustrated by the 9 x a diagrams of Figs. 15-17 for 8 — and p varying from p — 0.2 to p — 0.4, 
where there is at most only one tricritical point. This happens for all 8 < 0, and the tricritical point exists 
only for p < p mc , where p mc — > as 8 — > —00. Among the cases shown in the figures, the most interesting 
diagram corresponds to p = 0.3, where there are two points P C 2 and P' c2 at which 9t, 9 / and 9 C lines intersect. 

Figures 18 and 19 show a x 7 isotherms and 7 x 9 phase diagrams for 8 — with p = 0.2 and p = 0.3, 
respectively. In the first case, Fig. 18(b) shows the 7 x 9 diagram for a — —1.5 > a re , and three first-order 
lines are observed, one of which corresponds to 7 = 0. The two remaining symmetrical lines are tangent 
to each other for a = a re , and this gives rise to reentrant behaviour of the spontaneous magnetization for 
a < a re , as previously discussed. For a = a C 2 ~ —1.6549, as shown in Figs. 18(c) and (d), the critical 
temperature 9 C of the first-order 7 = line coincides with Q'p the temperature at which the symmetrical 
7^0 lines reappear. In the case p = 0.3, Figs. 19(b) and (d) exemplify the behaviour of the system around 
point P' c2 (Fig. 16), which marks the reappearance of a critical (second-order) point in the 7 = first-order 
line. Such critical point is absent for a' c2 < a < a C 2- 

5 Conclusions 

In this paper the one-dimensional Ising model with uniform infinite-range interactions and random nearest- 
neighbor interactions in an external field was studied. The random bond Ki connecting the sites i and i + 1 
was considered to satisfy the bimodal annealed probability distribution 

p(Kj) =p8(Kj - J A ) + (1 -p)8(nj - Jb). 

Explicit results were presented for the cases where there is competition between short- and long-range in- 
teractions. In order to avoid complete frustration of the system long-range interaction was assumed to be 
ferromagnetic. Competition was then produced by assuming Jb to be antiferromagnetic and letting J a be of 
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arbitrary character. It was shown that at T = the system can present antiferromagnetic, ferromagnetic and 
"ferrimagnetic" ordering. The spontaneous magnetization was shown to present unusual behaviour such as 
first-order ferri-ferromagnetic transitions, positive thermal derivative and reentrant behaviour. For certain 
distributions in which < 5 < 0.1649, where 5 is the ratio between the possible values of antiferromagnetic 
nearest-neighbor interactions, the system can present up to three tricritical points as the intensity of the 
infinite-range interactions is varied. Many topologically different field-temperature phase diagrams are pos- 
sible, and these can present up to four critical points, as opposed to the pure model which presents at most 
two critical points. 
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Figure captions 



Figure 1. a x 7(7 = /i/ | Js|) isotherms (a,c) and 7 x 8(9 = ksT / \Jb)\ phase diagrams (b,d) for 5 = | 
(<5 = Ja/Jb) and p = 0.5, with a = — 0.5(a = I/Jb) (a,b) and a = —0.6 (c,d). Dotted curves in (a) and (c) 
correspond to phase coexistence boundaries, while full lines in (b) and (d) indicate first order transitions. 

Figure 2. 9 x a diagram for 6 = § and p = 0.5. The dashed line shows the first order phase transition 
temperature (9 t ), which intersects the upper solid line corresponding to second order transition temperatures 
(9 C for h = and 9 hc for h ^ 0) at the tricritical point P tr . Temperature 8 hcri at which first order lines 
intersect in 7 x 9 diagrams is shown in the dash-dotted curve. Point Phcr2 marks the intersection of that 
curve with another second-order curve. In this and all the following 9 x a diagrams labels "para" and "ferro" 
relate to regions in the diagram where at zero field the system is at paramagnetic and ferromagnetic phases, 
respectively. 

Figure 3. The same as in the previous figure for S = | and p — 0.6. Solid and dashed curves represent second 
and first-order paramagnetic transition temperatures, respectively. The dotted curve shows the first-order 
ferri- ferromagnetic transition temperature 9f. At point P cr 3, where Ohcrz, Of and two 0t curves intersect, 
paramagnetic, ferrimagnetic and ferromagnetic phases coexist. The label "ferri" relates to the region where, 
at zero field, the system is at the ferrimagnetic phase. 

Figure 4. The same as in the previous figure for 6 = | and p = 0.7. Differently from the previous diagram, 
there is no three-phase coexistence point P cr3 nor point Phcr3, since the Of curve ends at P cri reaching the 
9' hc curve before reaching the t curve. 

Figure 5. Renormalized free energy / (/ = f a / \ Jb\) as a function of the spontaneous magnetization cr for 
various temperatures in the case 8 = |, p = 0.7 and a = —1.10 (a) and a = —1.13 (b). 

Figure 6. Spontaneous magnetization ao as a function of the renormalized temperature 8 for 5 = \,p = 0.6 
and various values of a. 

Figure 7. The same as in the previous figure for 5 — | and p = 0.7. Notice that for a = a cr there is a 
point in which ao is a continuous function of the temperature but the derivative diverges. 

Figure 8. 7 x 9 phase diagrams for 5 = |, p = 0.6 and values of a ranging from a = — 1 (a) to a = —1.21 
(f)- 

Figure 9. a x 7 isotherms (a,c) and j x 9 phase diagrams (b,d) for 5 = | and p = 0.7, with a = — 1 (a,b) 
and a = a cr ~ —1.1018 (c,d). 

Figure 10. 9 x a diagrams at zero field for 5 = \ and p = 0.4 (a), p = 0.48 (b), p = 0.5 (c) and p = 0.6 (d). 
Notice the reentrant behaviour of the first-order ferri- ferromagnetic transition lines, evident in (b) and (c). 

Figure 11. Spontaneous magnetization cr as a function of the renormalized temperature 9 for S — |, 
p = 0.5 and various values of a. 

Figure 12. (7x7 isotherms (a,c) and j x 9 phase diagrams (b,d) for 5 = g and p = 0.5, with a — a re ~ 
-1.2710 (a,b) and a = -1.29 (c,d). 
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Figure 13. 9 x a diagram for 5 = 0.16 and p — 0.453. Notice the existence of three tricritical points 
(P tr , P'tr an d P"r) an d of another zero-field second-order paramagnetic transition for a' t ' r < a < a' tr . For 
a' tr < a < a tr there is a third 9h c line just above the 9 t line, but omitted from the figure for clarity. 

Figure 14. 9 x a diagram for 8 = 0.15 and p — 0.42. The first-order ferri-ferromagnetic transition line 
(9f) intersects the second-order ferri-paramagnetic transition line (0' c ) and the first-order ferro-paramagnetic 
transition line (9 t ) at point P C 2, which is not a tricritical point. 

Figure 15. 8 x a diagram for 5 = and p = 0.2. 

Figure 16. 9 x a diagram for 5 — and p = 0.3. The dotted line in the inset corresponds to fist-order 
ferri-ferromagnetic transition temperatures (Of). 

Figure 17. 9 x a diagram for S = and p = 0.4. 

Figure 18. a x 7 isotherms (a,c) and 7x6* phase diagrams (b,d) for 8 — and p = 0.2, with a = —1.5 
(a,b) and a — a C 2 — —1.6549 (c,d). 

Figure 19. a x 7 isotherms (a,c) and 7 x 9 phase diagrams (b,d) for 8 = and p = 0.2, with a = —1.65 
(a,b) and a = —1.71 (c,d). 
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